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MANUAL PROCEDURE FOR DETERMINING POSITION IN 


SPACE FROM ONBOARD OPTICAL MEASUREMENTS 

By Harold A. Hamer 
Langley Research Center 


SUMMARY 


A method is developed for manually determining a position fix in space. 
Basically, the method employs four nonsimultaneous angular measurements of 
celestial bodies, along with precalculated data for a nominal trajectory, to 
determine vehicle position. The method is simple and accurate, with all onboard 
calculations required for one position fix being adequately performed with a 
slide rule in less than 15 minutes. 

The various equations required by the navigator, as well as those essen- 
tial for the precalculations, are outlined in detail. Although the method is 
developed and applied to Earth-Moon space, it could serve equally as well for 
interplanetary space. The important details of the method discussed are the 
linearity characteristics, the effect of star selection on accuracy, and the 
procedure to convert the nonsimultaneous measurements to a common time. 

The method is shown to have a relatively high degree of linearity; there- 
fore, it does not require the actual injection time to be close to a nominal 
injection time for good accuracy in determining position. On the basis of an 
error analysis of the method, the root -sum- square position error over most of 
the Earth-Moon distance is shown to be about 35 kilometers for a standard- 
deviation error of 10 arc-seconds in the angular measurements. Also, in rela- 
tion to the accuracy characteristics, almost any set of stars selected for the 
measurements will give good results, providing that the star used in the range 
determination lies within ±30° of the Earth-Moon- vehicle plane. 


INTRODUCTION 


In manned space missions the navigation and guidance will normally be 
accomplished by automatic or semiautomatic procedures. (For example, see 
refs. 1, 2, and 3 -) The inclusion of manual procedures (for example, ref. k) 
for these operations would be a desirable and, in some cases, a necessary fea- 
ture. With a manual system the astronaut could supposedly perform the essential 
parts of the midcourse navigation and guidance in the event of a malfunction 
in the space computer or power supply. Also, the manual system could provide a 
check on some of the results produced by the more complex automatic systems. 



The main requirements for any manual system would be simplicity, good accuracy, 
and fast calculation procedures. 

As a part of the overall navigation and guidance problem, a method is 
developed in this report which meets the requirements for manual position fixing. 
Basically, the method employs four nonsimultaneous angular measurements (taken 
by the same type of optical-measuring instrument) to determine range from a body 
center and then to determine the vehicle-position coordinates. The method is 
based on the astronaut* s possession of nominal-trajectory information and on 
certain precalculated data retained either in chart or table form. The astro- 
naut can adequately perform all necessary calculations with a slide rule and a 
scratch pad. 

In the present report the method is applied to an Earth-Moon trajectory, 
although it could serve equally as well for interplanetary trajectories. Impor- 
tant details discussed are the linearity characteristics of the method and the 
effect of star selection on accuracy. The procedure to convert the nonsimulta- 
neous measurements to a common time is also given. 

An error analysis is performed to show the accuracy characteristics of this 
method. The accuracy, of course, cannot be expected to approach that of auto- 
matic systems which use statistical reduction of a large number of various types 
of nonsimultaneous measurements. 


SYMBOLS 


A angle formed at vehicle by line of sight to Earth center and line of 

sight to Moon center 

AA,AB,AD,Z£ difference between actual values and nominal values 


B angle formed at Earth center by line to vehicle and line to Moon center 

b magnitude of vector perpendicular to plane containing r ev and h 

C angle formed at Moon center by line to vehicle and line to Earth center 

c constant (see eq. (ll)) 

D,E,F measurement quantities (see eqs. ( 15 )) 

h magnitude of vector perpendicular to instantaneous Earth-Moon- vehicle 

plane 


2, m,n direction cosine of line of sight to star with respect to X-, Y-, 
and Z-axis, respectively 

P,Q,R,S constants for any selected set of three stars (see eqs. (l7)\ 
r range or distance, (x^ + y 2 + z 2 ) ^ 
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Ar incremental range, r a -- r n 

rss root sum square (used as measure- of position-fixing accuracy), 

[("Ax ) 2 + (%r ) 2 + (“Aa) 2 ] 1 ^ 

R b body radius 

T time from injection 

t time 

At increment in time between common time (time for 11 common" or representa- 

tive value measurement) and time at which measurement is taken 

X, Y, Z vehicle- centered rectangular right-hand axis system in which X-axis 

is in the direction of Aries, XY-plane is parallel to Earth equato- 
rial plane, and Z-axis is in the direction of north celestial pole 

X r ,Y r , Z r rotating rectangular right-hand axis system in which X r -axis lies along 
Earth-Moon line, X r Y r -plane is in Earth-Moon plane, and Z r -axis is 
in northerly direction 

x,y, z position coordinates in rectangular right-hand axis system 

x r ,y r , z p position coordinates in rotating rectangular right-hand axis system 


x,y, z velocity coordinates in rectangular right-hand axis system 

Ax,£y,Az position off nominal trajectory in direction of x-, y-, and z-axis, 
respectively 

a one-half of angular diameter of body as viewed from vehicle 

5 angle formed at vehicle by line to star and its projection in the 

instantaneous Earth-Moon- vehicle plane 

0 angle formed at vehicle by line of sight to star and line of sight 

to body center 

0 angle formed at vehicle by line to Earth and projection of line to 

star in the instantaneous Earth-Moon- vehicle plane 

a standard-deviation value 


Subscripts: 

1, 2, 3 star 1, star 2, and star 3 
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specific stars 
actual (measured) value 


I, III, IV 

a 

A angle formed at vehicle by line of sight to Earth center and line of 

sight to Moon center 

m 

em distance between Earth center and Moon center 

ev distance or range from vehicle to Earth center 

mv distance or range from vehicle to Moon center 

n nominal value 

r range 

At incremental range, r a - r n 

unc uncorrected 

x,y, z component in direction of X-, Y-, and Z-axis, respectively 

Ax,Ay,Az position off nominal trajectory in direction of x-, y-, and z-axis, 
respectively 

0 angle formed at vehicle by line of sight to star and line of sight 

to body center 

Notation: 


CT 1 

O 


inverse of square matrix 
column matrix 


absolute value 

Bar over a symbol indicates a vector. 


- BASIC METHOD 


There are many combinations of optical measurements which can be used for 
position fixing. For example, reference 5 shows 13 different combinations of 
angular measurements which incorporate sightings on the Earth, Moon, Sun, or 
stars to yield solutions of vehicle position in Earth-Moon space. The basic 
requirements for a manual position-fixing procedure are that is must be simple 
enough to permit rapid calculations (several minutes) and still give accurate 



results. Of the various combinations of measurements, the only one that would 
seem applicable to this case is one that uses star-to-body- center angular meas- 
urements, together with a range measurement, as illustrated in figure 1. 




where Z, m, and n are the direction cosines of the respective stars. It is 
readily seen that the set of equations (l) can be, for the most part, precal- 
culated since the values of Z, m, and n are known in advance (by selecting 
certain stars). The solution of the set of equations (l), however, would 
involve calculating numbers up to six significant figures for Earth-Moon space. 
This problem can be circumvented by employing a nominal (or reference) trajec- 
tory for which vehicle position and various angles are known for any time along 
the trajectory. In this approach, the incremental values Ax, and Az 

(distances off the nominal) could then be easily and accurately calculated from 
the differences between the measured values of the angles and the values per- 
taining to the nominal trajectory. The equations in this form and some repre- 
sentative accuracies are shown in subsequent sections. 


5 


As shown by the set of equations (l), the vehicle position can be deter- 
mined by four simultaneous measurements, those of r, 0q, 0 2 , and 0^. In 

any practical case, simultaneous measurements are not possible, but a simple 
procedure can be used to convert nonsimultaneous measurements to a common time, 
as is subsequently shown. The range r in the set of equations (l) could be 
obtained by measuring the angular diameter of the body, that is, by replacing 
the quantity r with Rq, y4in a in which R^ is a known value of the body 
radius. However, this type of measurement is relatively inaccurate at large 
distances from a body and in Earth-Moon space can lead to position errors on 
the order of several hundred kilometers, as shown in reference 6. Therefore 
some other method must be used to obtain a measure of the range. For lunar 
trajectories, a triangulation procedure based on the Earth-Moon- vehicle triangle 
seems to be the most simple and accurate. Here again, this procedure is to be 
used in conjunction with a nominal trajectory for which the range is known for 
any given time along the trajectory. The details are presented in the following 
section. 


DERIVATION OF EQUATIONS 


Determination of Range 

In the present method for determining vehicle position along a lunar tra- 
jectory, a necessary first step is the determination of range to the center of 
the Earth. 

In figure 2 the angle ©q, which can vary between 0° and l80°, lies in the 
instantaneous Earth-Moon- vehicle plane. If figure 2 represents the vehicle 
position on a nominal trajectory at a given time, then the incremental range 

To star 1 
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Ar ev (defined as the actual range minus the nominal range) at that particular 
time is obtained as follows: 




ev sin A 


em sin C = . GIQ " - sin(A + B) 

sm A 


or 


r ev = r em cos B + r em sin B cot A 


( 2 ) 


As shown in equation (2), the variables which define r ev at a given time 
are the angles A and B; therefore: 


^ r ev . Sr ev 

dr™. - — — dA + — dB 

e SA 


Sb 


(3) 


From equation (2). 


Sr, 


ev 


Sa 


= -r em sin B csc^A 


but 


so that 


where 


sin B = sin A 
r em 


S^ev _ r mv 
5A sin A 


cos A = 


L ev 


+ r n 


mv 


- r 


em 


2 r r 
^ev^v 


(4) 

(5) 

(6) 


Furthermore, from equation (2), 


dr 


ev 


dB 


r e m(cos B cot A - sin B) 


(7) 


where the angle B is given by equation (4) and the angle A is given by 
equation ( 5 ). The values r 0V , r fflV , and r em are known values for the nominal 

trajectory. 
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The variables A and B in equation ( 3 ) refer to angular measurements 
which must be made onboard the vehicle. The angle B cannot be measured from 
the vehicle; therefore, an angle which is related to B and which can be meas- 
ured from the vehicle must be substituted. As shown by figure 2, the angle 0j_ 

measured between a star and the Earth center can be used in the following manner. 

For the case in which the star is in the instantaneous Earth-Moon- vehicle 
plane, d©i/dB = ±1, where the sign is determined by the relative directions of 
the Earth and the star. Equation ( 3 ), therefore, can be written as 

^4^4^ (8 > 

where d0q represents the change in the star-to-Earth angle and c is either 
±1. The equation for determining c is given in appendix A. If the star is 
not in the instantaneous plane, as shown in figure 2, d0p can be determined 
from the right spherical trigonometric relationship, 


cos 0j_ = cos 0^ cos 


or 


= arc cos 


/cos 0p\ 
1 cos Sjj 


Thus, 


d0p = 


sin 0 


cos^b^ - cos 2 0^ 


— — d0 -| 

1/2 1 


or 


1 - cos^B- 


1/2 


d0q = 


p p 

icos - COS 0J_ 


d0- 


( 9 ) 


Therefore, in terms of two angles which can be measured from the vehicle (Earth 
center to Moon center and star-to-Earth center), equation (8) can be written 


dr ev 8r evj 

dr ev = — r— — dA + c eV| 


cos 




1/2 


dA 


dB 


cos2s-^ - cos20jj 


d9- 


( 10 ) 
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or, 

and 


is linear with changes in A 


for the region in which a change in r ev 

B, 


^ r ev 

Ar ev = AA. + c 


5r ev/ 1 “ cos ^l 


,V 2 


8B Icos^Si “ cos^Bji 


AS- 


CII) 


. where 


AA. - A a - A n 


and 


- ®1, a “ 


%n 


(o° < Sjl < l 80 °) 


The incremental range Ar eV j then, is determined by the two 
and 0^ a and by the precalculated values of A n , 0^ n . 


measurements 
Sr e y Sr e y 
cA ' dB 


Aa 

and 


f 1 - cos^On \ ^ 

c — j . Equations for calculating values of 0, 5, and c 

\cos 2 6p - cos^G 2 _ ) 

which pertain to a nominal trajectory are given in appendix A, 


Linearity characteristics . - Investigation of the variations of Ar ev with 
AA and AB for a typical lunar trajectory showed that at a given time the 
partials br ev fbA and 5r ev /cffi can be considered to be essentially linear 

over a wide region, as indicated in figure 3* The nominal trajectory used 
( 70.6 hours from injection to perilune) was selected from reference 7 and is 
shown in figure 4. Table I presents pertinent position and velocity informa- 
tion of the trajectory, as well as values of the angles A and B. In fig- 
ure 3 the curve for AB was determined from the relations (see fig. 2) 


sin C 



sin A 


(12) 


and 


B = l80° - (A + C) (13) 

where A is held constant at the nominal value and is given by equation (6). 
The curve for AA was determined from the relation 

A r em . 
sm A = sm B 

r mv 
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4 x 10 U 



Figure 3.- Linearity characteristics of 
^ r ev ^ r ev 

and for a nominal dis- 

hA dB 

tance of l6l, 9^1 kilometers from cen- 
ter of Earth to vehicle. 


where 



Figure 4 .- Nominal trajectory. in rotating 
coordinate system; T denotes hours 
from injection. 


r mv 



-em 


2r r cos 
ev em 


B 


1/2 


and B is held constant at the nominal value and is determined by equations (12) 
and (13). The curves shown in figure 3 are for a nominal time from injection 
of 16.125 hours (r ev « 162,000 kilometers). At this time a 2- hour difference 

at injection from the nominal injection time would represent a difference in 
range from the Earth of about 13*000 kilometers. The data in figure 3 show 
that this difference can be considered to be within the linear region. Several 
nominal or reference trajectories, therefore, could suffice for range determina- 
tion for a given daily launch window. 

The linearity characteristics of dr ev /c)A and br eY fc)B illustrated in 

figure 3 are representative of those over most of the Earth-Moon distance. 
Actually the linear approximation for Sr ev /dA improves as the Moon is 
approached, as shown by the derivative of equation ( 5 ) which is 


6 r ev 
6A 2 


r mv 


cos A 
sin^A 
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TABLE I.- NOMINAL TRAJECTORY PARAMETERS 


[vehicle position and velocity coordinates are vehicle centeredj 





Time from injection, hours 



Variable 









1 

2 

3 

4 

5.5 

IO.5 

16.125 

Xgyj km ...... 

17, 000.141 

19,873-013 

20, 772.521 

20,874.339 

20,284.239 

15 , 797.195 

9,255-8512 

yg'yj km ...... 

13,826.706 

28,971.665 

41 , 465-015 

52,1+16.910 

66,96 5.325 

106,257.19 

141 , 266.32 

z gyj km •••••• 

8,591.5583 

17 , 11 + 1.773 

24 , 114.484 

30,189.287 

38,219.273 

59,747.541 

78,788.545 

r } km 

25,537.152 

39,091.360 

52,271.918 

63,989.558 

79 , 727.772 

122,922.38 

162,016.91 

km ...... 

-189,769.15 

- 183 , 757-26 

-179,697.19 

-176,1+13.57 

- 172 , 192.03 

- 160 , 316.94 

- 147 , 900.10 

7mv> ...... 

-254,284.19 

-241, 070.ll 

- 230 , 478.64 

-22^,399.1+3 

-209, 604-99 

-179, 008.56 

-152,872.43 

Z^py, km ...... 

-129,595.21+ 

-122,195.85 

-• 116 , 359-00 

- 111 , 404.98 

-105,027.68 

-88,757.164 

- 75 , 157.655 



31+2,735.63 

326,823.12 

3 li+, 56 i +.62 

304 , 221.18 

290,886.84 

256,170.68 

225 , 594.98 

r em' ^ 

365, 691+. 67 

365 , 790.69 

365,888.32 

365 , 987.50 

366 , 139.23 

366,669.32 

367,307.81 

x ev , km/sec .... 

1.3791860 

O.42971609 

0.11338018 

- 0.041338963 

-0.16327614 

- 0.30064959 

- 0.33717387 

y ev , km/sec . . . . 

1+. 7862262 

3.7612012 

3.2252961 

2.8804782 

2.5318916 

1.9188737 

1.5703809 

z gyj km/ sec .... 

2.7279185 

2.1078791 

1.7935755 

1.5940144 

1.3941551 

1.0468713 

0.85165680 

x^., km/sec .... 

a 2.2l+ 

a l .67 ■ 

a 0 . 973 

a 0 . 84 l 

a 0.746 

a 0.626 

a o. 6 i 7 

y mv , W sec • - • • 

a l+.l+0 

a 3-67 

a 2-70 

a 2.38 

a 2.08 1 

a l.50 

a l.l5 

Zj^^ km/ sec .... 

a 2.38 

a 1.82 

a 1.46 

a 1.28 

a 1.12 

a 0.788 

a 0.6o4 

A, deg ...... 

i 66 ° 3055 " 

175°650” 

168 ° 725 " 

163°1262" 

158°958" 

I 48 °l 88 l" 

142 ° 774 " 

B, deg 

12°1125" 

4°1099» 

10 ° 449 ” 

13°2801" 

17 ° 392 " 

2i°l4o4" 

22 ° 379 " 



23.125 

31.125 

45.125 

60 

65 

68.5 

69.5 

km *•«••• 

618.92375 

- 9 , 233.1921 

-25,632.839 

-40,898.372 

.44,885.478 

- 46 , 341.619 

- 45 , 992.488 

y ev , km ...... 

177,258.93 

211,785.18 

261,074.69 

303 , 197.37 

315 , 953.95 

324,980.53 

327,676.34 

Zgyj kin •••••• 

98,256.097 

116 , 836.74 

143,211.87 

165,528.05 

172,234.22 

176,961.26 

178,384.29 

r0Y, km ...... 

202, 670.61 

242, 051.73 

298,868.83 

347,851.77 

362,637.86 

372,927.85 

375,909.62 

km •••«•« 

- 132 , 226.49 

-113,475.61 

-78,287.308 

-37,619.423 

-22,760.254 

- 11 , 053-984 

-6,951.3123 

ymv, km ...... 

- 126 , 541.42 

- 101 , 136.23 

- 62 , 762.242 

-25,101.552 

-12,194.307 

- 2 , 572.1041 

-367.61865 

Zjpy, km ...... 

-61,739.027 

-49,059.589 

-30,346.349 

-12,316.612 

- 6 , 164.1073 

- 1 , 558.2646 

-129.66503 

r mv' ^ 

193,153.52 

159,725.0 6 

104,827.97 

46,872.249 

26 , 546.686 

11,455.759 

6,962.2336 

r em-» ^ 

368,158.94 

369,198.98 

371,160.79 

373 , 389.42 

374 , 161.41 

374,706.77 

374 , 863.24 

x ey , km/sec .... 

-0.34476699 

-0.33782194 

-0.31034989 

a -o. 252 

a - 0.181 

a 0. 00410 

a o . 243 

y ev > W s ec . . . . 

1.3058873 

1.1039615 

0.87009046 

a 0.735 

a 0.704 

a 0.74o 

a 0 . 744 

z ey , km/sec .... 

0.70449040 

0.59265123 

0.46342536 

a 0.382 

a 0.369 

a 0.389 

a 0 . 398 

Xjjjry, km/ sec .... 

a o.632 

a 0.666 

a 0.725 : 

0.79926153 

0.86291157 

1.0488207 

1.2677780 

y mv , W sec • • • • 

a 0.964 

a 0.830 

a 0.720 

0.70492821 

0 . 73557 H 9 

0.80363455 

0.82276785 

Vv, W sec • • • • 

a 0.469 

a 0.405 

a 0.347 

0.33600494 

0.35090673 

0.38753199 

0.40476791 

A, deg 

136 ° 3204 " 

132°l683" 

126°2528 M 

119 ° 2995 M 

H 3 ° 3074 h 

98 ° 234 " 

8o°2992 n 

B, deg ...... 

21° 260" 

18 ° 2197 ” 

13 ° 316 ” 

6 ° 907 M 

3 ° 2592 ” 

1°2646" 

l°l 84 " 


approximate value is given inasmuch as this value was not calculated in the trajectory program used In the 
present analysis. 
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As noted from table I, the angle A approaches 90° as the Moon is approached. 
From the derivative of equation ( 7 ), 


= -r (cot A sin B + cos B) 

Bb 2 

it can be determined (by use of the values of A and B from table I) that 
the linear approximation for dr ev /cffi will also improve as the Moon is 

approached . 

Effect of proximity to Earth-Moon line .- The values of Sr ev /dA and 
Sr ev ^/cffi which are precalculated and supplied to the astronaut in the form of 

charts or tables are shown in figure 5 for vehicle positions along the ascent 
portion of the nominal trajectory. The partials dr mv /dA and c)r mv /c)C 

required to determine range from the Moon are also shown. These values will 
approach those shown for Sr ev /dA and Sr ev /dB at distances near the Earth. 

It is to be noted that for a position fix near the Moon higher accuracy is 
obtained by measuring range from the Earth rather than from the Moon, inasmuch 
as the geometry gives lower values for Sr ev /cA and dr ev /c)B. From equa- 
tion (ll) it is seen that the error in determining the distance off the nominal 
Ar is proportional to the magnitudes of the partials. 



Figure 5.- Variation over Earth-Moon distance in partials required for range determination. 
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In figure 5 the large values shown for 5r ev /dA and dr 0V /dB near the 

Earth occur because of the close proximity of the vehicle to the Earth-Moon 
line (i.e., the angle A approaches l8o°). The relative position of the vehi- 
cle with respect to the instantaneous Earth-Moon line is shown in figure 4. 
Thus, range or position-fixing determination by use of equation (ll) would be 
inaccurate for about the first 5 hours (or 75,000 kilometers) from injection. 
These results would be typical for any Earth-Moon trajectory for any method 
which is based on Earth-Moon measurements. The accuracy for the first 5 hours 
could be improved considerably by determining range from the angular-diameter 
measurement of the Earth, as shown in figure 6. "The dashed curve represents 
the error in range for an angular-diameter-measurement a error of 10 arc- 
seconds. 

For times along the trajectory except for the first 5 hours, the error in 
determining range by equation (ll) is relatively small, as shown by the solid 
curve in figure 6. This curve represents the error in the incremental range 
Ar ev along the ascent portion of the nominal trajectory where the angular 

measurements A and 0j_ are considered to have random uncorrelated errors, 
each with a standard deviation a of 10 arc- seconds such that 



Nominal distance from Earth center, r ev > km 


Figure 6.- Variation of error in range determination over Earth- Moon distance; a of angular- 

measurement error is 10 arc-seconds. 
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In figure 6, the star which is used to measure 0p is assumed to lie in (or 
near) the instantaneous Earth-Moon- vehicle plane so that the value of the term 


1 - cos^Gp 

p p 

icos^bp - cos ©ij 


1/2 


in equation (ll) is 1 and 


Sr 


ev 


SB 


= ± Ihc. 

Sg i 



Figure 7 .- Effect of star position on accuracy of range, determination; 

vV2 




^ev = 


hk 


■ Ak + c 




2 

cos 0 


cffi \cos^5 - cos^0 


£6. 


Effect of st ar selected for meas urement.- The value of the quantity 

*V 1 2 “" 


can vary from 1 to °o, depending on the nominal trajectory 


1 - cos^0p 
ycos^Sq - cos^Gp^ 

and the star selected for the measurement of 0p. The effect of star position 


on the error in Ar is shown in figure 7 . The error in Ar is a minimum when 

vl/2 

a 


the quantity 


1 - cos^G- 


cos 8l - cos 2 0p 


is a minimum. This quantity has its minimum 


with respect to 0 at 0 = 90 ° and rises rapidly as the extreme values, 0=5 
and 0 = l80° - 5, are approached. The extreme values of 0 occur when the 
star position is in the plane containing the vehicle-Earth line and the vertical 
to the Earth-Moon- vehicle plane. The desired stars, therefore, are those that 
lie near the Earth-Moon- vehicle plane (5 within about ±30°) with projections 
in that plane away from the- vehicle-Earth line. In figure 8 values of the 




(a) Angle between star and Earth center. 




Nominal distance from Earth center, r gv , km 

(b) Angle between star and Its projection 
in instantaneous Earth- Moon- vehicle 
plane . 


Figure 8.- Variation of 8 and 5 over the Earth-Moon distance for nominal trajectory. 


angles 0 and 6 for several stars are shown along the ascent portion of the 
nominal trajectory. 

As indicated in figure 8(b), large changes can occur in 5 as the vehicle 
crosses "under" the Earth-Moon line at r ev ~ k-0 , 000 kilometers. In this region 
the Earth-Moon- vehicle plane rotates about the Earth-Moon line through an angle 
of approximately l80°. For the star Regulus, which forms a star-to-body angle 
0 of about 90 °, the value of 6 would go. through a range of angles close 
to 180 °. 


Determination of Vehicle Position 

As previously stated, equations (l), in the form given, are not convenient 
for manual, calculation of vehicle position. The positions off the nominal 
trajectory, however, can be easily calculated from the following set of 
equations: 


x a - 

*n 


Z 1 

m i 

n l 

-1 

r a 

COS 

H 

CD 

a “ r n 

COS 

•l.n'l 

y a - 

y n 

) = 

z 2 

m 2 

n 2 

< 

r a 

cos 

0 2, 

a ” r n 

cos 

0J 

CD 

z a 

V. 

z n 

J 


b 

m 3 

n 3 


r a 

cos 

e 3. 

a - r n 

cos 

e 3,n 

J 


(14) 


where the subscript a represents an actual (measured) value and the sub- 
script n represents a nominal value. After substituting 

x a - x n 

y a - y n ) 

z a " z n 
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and 


^r a cos Q 1)& - r n cos 0^^ 

“ < r a cos 0 2 ^ a - r n cos 0 2#n V 
Ja cos 0^ a - r n cos 0^ 
a solution of the set of the three simultaneous equations (l4) yields: 

Az = 

Ay = 

Ax = 

where 

P = ^m^Zp “ m 2^l) 

Q = ( n 1^2 " n 2^l) ^ 

R = ( m l Z 3 “ m 3 2 l) 
s = (n^ - n^) 



I 

F 

V. J 


(15) 


(16) 


(17) 


The parenthetical expressions in equations (l6) and (17) are constants for any 
three given stars and can be precalculated. Various solutions of the set of 
equations (l4) can be obtained, but all would be similar to those given by 
equations (l6). 


Therefore, equations (l6), or some similar to these, are the equations that 
an astronaut could use to calculate vehicle position. The star-to-body angle 
measurement utilized in equation (ll) is used, of course, as one of the star 
measurements in the set of equations (14). As shown by equations (l6), a posi- 
tion fix would require, first, a calculation of the quantities D, E, and F 
from the four measurements r a , 0^ a , e 2, a-> an ^ a- 7 secon( l.> nine slide rule 

multiplications, and third, several additions and subtractions on a scratch pad. 
Inspection of the expressions for D, E, and F (eq. (15)), however, indicates 
calculations involving a large number of significant figures which cannot be 
accurately performed on a slide rule. A procedure similar to that previously 
shown for determining Ar can be employed by making use of precalculated values 
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of D, E, and F for a nominal trajectory, as described in the following 
section. 



Determination of the meas urement quantities D, E, and F .- The quantity 
D is illustrated in figure 9 . By inspection of figure 9, the derivative of D 
is seen to be: 


Sd 

dD = - dr, 


Sr, 


ev 


Sd 

+ aeT 01 


or, for the region in which a change in D is linear with changes in r ev 
and Op, 


Sd 

AD — — Ar ev 
or e v 


+ 


Se 1 


A9p 


Since D = 0 when r a cos 0q^ a = r n cos 0~^ n j then ZD = D, so that 


D 


5d 

Sr ev 


Ar e v 



(18) 


From the expression 


D = r* cos 0 


l,a 


r n cos ei,n 


it is seen that 


SD 


Sr 


ev 


- cos 0-j_ 


(19) 


* 
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and 




-r ev sin 6 -l 


( 20 ) 


■where r ev and 0-^ are known values for the nominal trajectory. Equations 

similar to equations (l8), (19), and ( 20 ) are obtained for E and F by 
replacing 0 1 with 0 2 and 9 y respectively. 

Values of the partials dD/8r ev and dD/a0 are shown in figure 10 for the 

example nominal trajectory. The data correspond to the star Regulus, which lies 
near the instantaneous Earth-Moon- vehicle plane throughout most of the trajec- 
tory. (See fig. 8(b).) The astronaut, then, would be supplied with informa- 
tion, such as given in figure 10 , concerning three preselected stars. This 
type of information, along with the values for Ar ev , an< ^ 

would allow him to calculate readily the quantities D, E, and F by equa- 
tions similar to equation (l8). The value of Ar ev is obtainable directly 

from equation (ll). 

Linear region for the quantities D, E, and F.- For equation (l8) to be 
exactly correct, the partials dD /8r ev and dD/dep for a given time must be 

constants. Inspection of equation (19) shows that dD/8r ev is a constant for 
all values of r ev . From equation ( 20 ), 


Hence, at a given time. 



dD/801 will 



-r 


ev 


cos 0j_ 


be constant for a value of 9 q = 90°. 
Actually, cffl/80p will be approximately 

constant over a large range of values of 
0 ^ (for example, 60 ° < 0 p < 120 °). Star 

selection should not be necessarily 
limited to this range of star-to-body 
angles, however, unless very large devia- 
tions (greater than ±2°) from the nomi- 
nal 0 are expected. This deviation of 
2° would correspond roughly to a 2-hour 
difference between the actual injection 
time and a nominal injection time. 


CORRECTION FOR NONSMJLTANEOUS 


Figure 10.- Variation over Earth- Moon dis- 
tance of and of f or B tar 

o0 drg-y- 

Regulus . 


MEASUREMENTS 


In the equations for determination 
of range (eq. (ll)) and position 
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(eqs. (16) ), all four measurements are 
assumed to have been made simultaneously. 
Since simultaneous measurements are 
usually not possible with the same 
instrument, errors would be introduced 
in the results unless the measurements 
were corrected or converted to some com- 
mon time. This time (for which the 
position fix is required) could be 
either before or after the actual time 
of the measurements. If the time cor- 
responding to one of the measurements is 
made the common time, only three meas- 
urements would have to be converted. 

The measurements can be converted 
to a common time (small increments of 
time, only) by use of precalculated data 
on the rate of change of angle with time 
along the nominal trajectory. Examples 
of such data are shown in figure 11 for 
several stars, as well as for the angu- 
lar measurement between the Earth and 
the Moon. (Table II identifies the 
various stars used in this report.) 

These data would be practically the same 
for any other trajectory reasonably near 
the nominal, so that they would apply 
for making corrections when off the 
nominal trajectory. It is of interest 
to note that the curve shown for dA/dt 
would be about the same for any trajec- 
tory to the Moon providing the flight 



Figure 11.- Rate of change of measure- 
ment angles with time along ascent 
portion of nominal trajectory, where 
A is angle between Earth and Moon 
centers and 0j, ©m, 0 jv are 

angles between different stars and 
Earth center. (See table II for 
identification of stars.) 


TABLE II.- STAB POSITIONS USED IN ANALYSIS 


Star 


m 


Number 

Name 

l 

n 

i 

Capella 

0.13862 

0.68133 

0.71873 

ii 

Sirius 

-.17870 

.94126 

- . 28652 

hi 

Procyon 

-A0895 

.90781 

.09295 

IV 

Regulus 

-.85903 

.46641 

. 21104 

V 

Rigil Kentaurus 

-•37937 

-.30996 

-.87178 

VI 

Vega 

.11866 

-.77099 

.62570 



time is in the neighborhood of JO hours. Equations for calculating dA/dt and 
de/dt along a nominal trajectory are given in appendix B. 


In converting a measurement to a common time, the corresponding value of 
d0/dt or dA/dt at the given range multiplied by the time increment would be 
sufficient (for converting over short time intervals) to calculate the incremen- 
tal angle to be added algebraically to the measured angle. The range would need 
to be known only approximately. For trajectories launched near nominal injec- 
tion time, the nominal range (at the time of measurements) would be adequate, 
whereas for trajectories not launched near nominal injection time the approxi- 
mate range could be determined according to time from actual injection time. 
Because the correction to the angle A would ordinarily be the largest (fig. 11 ) 
and also because the precalculations for dA/dt are more extensive and require 
more information on the nominal trajectory (see appendix B), it would be best 
to use the time at which A is measured as the common time. In this case. 


A a “ ( A a)unc 

a “ (®l,a) un c + at 

/x ae P ) 

^a v ^a/unc &t ^ 

d0 ■z 

0 3,a = ( e 3,a) unc + aT At 3 



Nominal distance from Earth center, r , km 

ev 


40 x 10^ 


Figure 12.- Accuracy of manual method over Earth-Moon distance; a of measure- 
ment error is 10 arc-seconds. (See table II for identification of stars.) 
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ERROR ANALYSIS OF METHOD 


The calculated results shown in figures 12 and 13 are typical of the accu- 
racy obtainable with the present method, providing the range is detennined by 
the angular measurements given in equation (ll). It can be seen that after the 
first 5 hours the root- sum- square position error is in the neighborhood of 
35 kilometers. The large position errors shown for the first 5 hours are due to 
large errors in range determination (fig. 6) when the vehicle is in close prox- 
imity to the Earth-Moon line. If, however, the range in this region is deter- 
mined by an angular- diameter measurement of the Earth, the root- sum- square 
position error could be reduced to 40 kilometers or less. 

In the analysis the errors in each of the four measurements were assumed 
to be uncorrelated. The results are for the case in which the star used in 
determining Ar lies in or near the instantaneous Earth- Moon- vehicle plane; 
that is, the range error shown in figure 6 would apply in this case. In the 
error analysis, linear perturbation theory was used to calculate values for the 
root sum square of the standard deviations of the component position errors. 

(See appendix C for the pertinent equations.) This root- sum- square value can 
be considered as a measure of the position error. 

The results in figure 12 are for standard deviations of measurement error 
of 10 arc-seconds in each of the measurements. It should be noted that doubling 
the measurement error will approximately double the position-determination 
error. The variation in position error with range is shown in figure 12 for 



Nominal distance from Earth center, r^> 


Figure 13 •- Comparison of manual method with two other methods of high accu- 
racy; a of angular- measurement error is 10 arc- seconds and of radar- 
distance-measurement error is 16 meters. 




cases in which two different sets of stars are used for the measurements. The 
accuracy for either set is about the same, and this fact would be true for any 
combination of stars selected, with three exceptions. These three exceptions 
are that (l) the star used for the determination of Ar (eq. (ll)) must lie 
within about ±30° of the instantaneous Earth -Moon- vehicle plane (fig. 7)> (2) 
the stars selected should be such that the rate of change of 0 with time along 
any part of the trajectory is comparatively large (ref. 8); for example, such as 
is indicated for the stars Regulus and Procyon in fig. 8(a); (3) the stars 
should not all be at approximately the same location. 

In figure 13 the manual method is shown to be about as accurate over most 
of the Earth-Moon distance as two other position-determination methods of good 
accuracy. The results shown for these two other methods were obtained from 
reference 6. The radar-measurement method consists of the same onboard star-to- 
body angular measurements as the manual method; however, the range is deter- 
mined by a highly accurate Earth-based radar measurement. The errors shown for 
this method are almost entirely due to the errors in the three angular measure- 
ments. In the five -angular- measurement method, five onboard angular measure- 
ments are used to determime vehicle position. The five angular measurements are 
the included angles between stars and body centers. More specifically, the 
measurements are the included angles between each of three stars and Earth (or 
Moon) center and between two of these stars and Moon (or Earth) center. Of all 
the optical methods which rely completely on onboard sightings and which do not 
use statistical procedures or information on a precalculated nominal trajectory, 
this five-measurement method is the most accurate. (See ref. 6.) However, 
because of the complexity of the equations, the onboard calculations for this 
method would require an automatic computer. 


SUMMARY 0E PROCEDURES FOR MANUAL METHOD 


No attempt is made in this report to systematize the present method. The 
required procedures are simple enough and any system used would be according to 
individual preference. The various required calculations (both precalculations 
and onboard calculations) are summarized in the following outline in the proper 
order, along with the pertinent equation numbers. 

Precalculations required for the nominal trajectory for various times from nom- 
inal injection time (assuming that nominal values of the quantities listed in 
table I are known): 


Angle A 

Angles 0 1> 9 2' and 0^ 

d9-j^ d0p d0^ 

Derivatives , — and — 

dt ' dt ' dt dt 

dr dr 

Partials and — — 

dA dB 


.... Equation (6) 
.... Equation (Al) 

Equations (Bl) and (B3) 

Equations ( 5 ) and ( 7 ) 
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1/2 


Quantity 


Partials 



by use of Equations 


5e 


and 


3b ? 


' So ’ 

' ev w "ev "ev W 1 
According to the selected set of stars, calculate the 
constants P, Q, R, and S (eqs. (17)), as well as the 
other constants appearing in equations (l6) 


Equations 


(A2) and (Ab) 
( 19 ) and (20) 


Onboard calculations required for a position fix (assuming that the four non- 
simultaneous measurements have been taken over a relatively short time 
interval) : 

Convert each measurement to a common time (small increments 


of time, only) by equations similar to Equations (21) 

Determine the values AA, AB^j £&2> an< ^ are the 

differences between the actual and nominal angles 

Calculate the value Ar ev Equation (ll) 

Determine the quantities D, E, and F from equations 

similar to Equation (l8) 

Calculate position Ax,Ay,Az relative to nominal 

trajectory Equations (l6) 

Calculate true position x,y, z by algebraically adding Ax, 

Ay, and Az to nominal trajectory values 


The units for the onboard calculations would be arc- seconds for the angles 
and kilometers for the range and position values. 


A trial calculation of a position fix which would apply to an Earth-Moon 
trajectory was performed in approximately 15 minutes. For this calculation it 
was assumed that the measurements had been recorded. With two persons working 
together, the time could be reduced. 


CONCLUDING REMARKS 


A simple method for manual determination of position of a space vehicle 
has been presented. The method as presented has been applied to a lunar trajec- 
tory; however, it would apply equally as well to an interplanetary trajectory. 

The method makes use of the knowledge of the distance between the Earth 
and the Moon (or the distance between a planet and the Sun) to obtain an accu- 
rate measure of position relative to a precalculated nominal or reference tra- 
jectory. The linearity characteristics of the method are such that good accu- 
racy can be obtained even though differences as large as 2 hours exist between 
actual and nominal injection time. In addition to the one type of instrument 
to measure angles, the only other requirements for the navigator are time- 
history information on the nominal trajectory, various charts or tables of 
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precalculated information, and several simple equations which can be solved 
adequately with a slide mile. It has been determined that after the measure- 
ments have been made the method would require no more than a total of 15 minutes 
to calculate a position fix. 

From results of an error analysis, the accuracy of the method after about 
the first 5 hours from the Earth was shown to be approximately constant over 
the remaining Earth-Moon distance. In this region the error in position deter- 
mination was shown to be in the neighborhood of 35 kilometers. For the times 
within the first 5 hours, close proximity of the vehicle to the Earth-Moon line 
causes large error in position determination unless range is determined by a 
measurement of the angular diameter of the Earth. Also, in relation to the 
accuracy characteristics, practically any set of three stars selected for the 
measurements would give good results providing that the star used in the range 
determination lies within about ±30° of the instantaneous Earth-Moon- vehicle 
plane . 

Under linear assumptions, a manual guidance procedure could be developed 
which makes use of the position measurements. The same number of calculations 
as for two position fixes could determine the velocity changes for performing 
a midcourse correction. The main limitation to such a procedure would be the 
distance the flight trajectory could be allowed to deviate from the nominal 
trajectory. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 20, 1964. 
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APPENDIX A 


EQUATIONS FOR CALCULATING 6, 0, AND c 

The angles 6 and 0 (fig. 1*0, as well as the value c, pertaining to a 
nominal trajectory can be precalculated for any given star with known direction 
cosines 7, m, and n by the following equations. 

To star 



Figure lk.- Illustration shoving the vectors b and h. 

The angle 0, which is the angle from the vehicle between the lines of 
sight to a star and a body center (Earth, for example), is obtained from the dot 
product of the position vectors of the Earth center and star 


cos 9 = ZX ev + my ev + . n .feY ( 0 ° 1 9 ^ l80°) <A1) 

r ev 

By defining a vector h perpendicular to the instantaneous Earth-Moon- vehicle 
plane, the angle 6, which is the angle between the line of sight to the star 
and its projection in this plane, is obtained from the dot product of the vec- 
tor h and the position vector of the star 


sin 6 


2h x + mhy + nh z 

_ ~ h ~ 


(A2) 


where the absolute value can be used such that 0° ^ 6 ^ 90 ° and where the 
values for h x , hy, h z , and h are obtained from the cross product of the 
position vectors of the Earth and Moon centers. Thus, 


25 



(A3) 


Ax = Yev^nv ~ z evymv 

A y — z ev x mv “ x ev%T / 
A z = x ev y mv - y eY x mv 


and 

/ p 2 2 

A = \h x + hy + h z 

The sign of" the quantity c in equation (ll) is determined by the posi- 
tion of the star with respect to the Earth-vehicle line r e v The, sign is 
positive (negative) if, as viewed from the direction of h, the projection of 
the line to the star in the instantaneous Earth-Moon- vehicle plane is to the 
right (left) of the Earth- vehicle line. As shown in figure 14, the relative 
directions of these two lines at a given time along a nominal trajectory can be 
determined by first defining a vector b which is perpendicular to the plane 
containing the vectors r e y and h. 

The vector b is obtained from the cross product of the vectors r e v 
and h such that the components of b are given by 

Ax = yevAz “ z ev^y 


1/2 


by — z ev h x - x ev h z 


A Z — Xg-yhy 


yevAx 


where the values for h x , h y , and h z are given by equations (A3). Next, the 
dot product of the position vector of the star and the vector b will yield 
the component of the star vector along the vector b which is given by the term 
2b x "t* mb v nb z 

. If this term is positive, the projection of the star vector 

b 

in the instantaneous Earth-Moon- vehicle plane is to the right of the Earth- 
vehicle line r ev so that the quantity c is given as 


Zb x + mby + rib z 
7b x + mby + nb z 


( Ah ) 
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APPENDIX B 


EQUATIONS FOB CALCULATING RATE OF CHANGE OF 
MEASUREMENT WITH TIME 


The rate of change of the measurement 0 (angle between star and Earth 
center) with time along a nominal trajectory is given by 

d0 50 ^ev , 50 ^ y ev , 50 ^ z ev 

— + + 


dt 


5X0-y dt 


5y ev dt 5z ev dt 


(Bl) 


where dx ev ^dt, dy ev /dt, and dz eV /dt are known values for the nominal trajec- 
tory and, as can be determined from equations given in reference 8, the partials 
for the angular measurement 0 are 


50 


c ev ( m y ev ' — ev 


+ nz„ 


)-*( 


y e v 2 + z ev 2 ) 




ev 


dx ( 


ev p 

r ey 


r ev 2 - ( Zx ' 


ev + ^ev + nz ev 


) 2 ; 


1/2 


50 


y ev (^ x ev + nz ev 


) - m ( 


2 

x ev + z ev ) 


dy 


ev 


„ 2 „ 2 
r ev r ev 


(Zx 


ev + ^ev + nz ev 


) 2 ] 


1/2 


(B2) 


5© 


Sz 


ev 


^ev 

2 


(zx ev + ^ev) " n ( x ev + ^ev ) 

1/2 


■ev 


_ r ev 2 - ( ! *e 


+ my PV + nz <= 


o 2 ] 


The rate of change of the measurement A (angle between Earth and Moon 
centers) with time along the nominal trajectory is given by 

dA _ 5A d^ev + 5A ^ev + 5A <3 ~ z ev + 5A ^^mv + 5 a ^mv 
dt 5x ev dt 5y ev dt 5z ev dt 5x mv dt 5y mv dt 


5A ^■ z nrv 


Szjuv dt 


(B3) 


where dx ev /dt, dy ev /dt, dz ev /dt, dx mv /dt, dy^/dt, and dz mv /dt are known 
values for the nominal trajectory. The partials of A with respect to the 
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Earth are determined from equations (B2) by replacing the direction cosines of 
the stars l, m , and n with the direction cosines of the Moon: 

5A _ x ev{^ev^mv + z e v z mv) “ x mv(yev^ + z ev^) 

Sxev ?r p 2 , 2I 1 / 2 

r ev jj^ev r mv ” \ x ev x my y ev^mv' + z ev z mv) J 
Sa _ y ev ( X ev x mv + z ev z mv) - ^mv ( x ev + z ev ) 

^ev of 2 2 / vol 1 / 2 

r ev [Jev r nrv ~ ( x ev x mv + ^ev^mv + z ev z mv) J 

5A _ z ev( x ev x mv + ^ev-^mv) “ z mv ( x ev + ^ev ) 

^ Z ev P -.1/2 

r ev [Jev r znv ~ ( x ev x mv + ^ev^mv + z ev z mv) J 

The partials of A with respect to the Moon are determined from equa- 
tions (B2) by replacing x ev , y ev , z ev , and r ev with y mv , z^, and 

r mv^ respectively y and by replacing the direction cosines of the stars l> m, 
and n with the direction cosines of the Earth: 

5A _ x mv(y e yymv + z e v z mv) ~ x ev(ymy 2 + z my 2 ) 


r ev r mv " ( x ev x mv + y 


ev^mv + z ev z mv) 


5A _ ymv( x ev x mv + z ev z mv) “ yev( x mv + z mv ) 

^mv 2 T 2 2 Al 1 / 2 

r mv jjev r mv "" ( X ev x mv + ^ + z ev z mv) J 

5A = + yeyymv) ~ z ev ( x my 2 + y^ 2 ) 

^ Z mv i“ “I 2 

r mv |^ev r mv “ ( X ey x niv + ^ + z ev z mv) J 
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APPENDIX C 


EQUATIONS FOR ERROR ANALYSIS OF METHOD 


For a linear perturbation procedure, which is used in this report to 
analyze position error, the system of equations (l) differentiated with respect 
to 9 ^ yields equations for determining the effect of errors in 6^ on the 

position coordinates x, y, and z: 


Z 1 



m l 



dz 

1 a? 


-r sin 0^1 


SB i 


+ m^ 


+ n P — = 
^9 c50 


Sx Sy c)z 

__ + __ + _ = ° 




(Cl) 


When linear perturbation is used, the partials are, of course, assumed to be 
constants. Data in reference 6 show that these partials are essentially con- 
stant over variations of 0 of ±1.5°- These changes in 0 would represent 
trajectories which are thousands of kilometers off the nominal trajectory. 


A solution of the simultaneous equations (Cl) yields: 

3z _ r sin e 1 (PZ 3 - RZ 2 ) 

&1 " QR - PS 

3y _ ~ l 2 r sln e i Q 3z 

30! 1 p 3e7 

3x _ mp 5y ng 3z 

ig 3 g x Zg Se-L 

where P, Q, R, and S are defined by equations (17). 

In a similar manner, the partials for 0g, 0 3 , and r can be determined 
as 


RZjjr sin 0g 
30^ = QR - PS 
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5y _ S 5z 

g r 5e 2 

_ ™i 5y n i 5z 

5b 2 5©2 l± 5©2 

5z _ Pill* sin 0J 
5© 3 OF - PS 

dy _ _ Q 5z 
503 P 5©3 

5x _ m i 5y n l 5z 

503 2i 503 l l 50 5 

cos 0 j_^RZ2 - P^ 3 ^ + ix^P c os 03 - R cos Bq') 

QR - PS 

cty m l _2 cos 9 1 - *1 cos 9 2 Q 5z 
5r P "" " ’ “ P 5r 

5x = cos 9 1 _ ^ 5y _ ^1 5z 
5r i 1 5r l 2 5r 

The standard deviations of the errors in the x-, y-, and z-directions due 
to errors in the measurements are, then. 



The root- sum-square (rss) of the position error due to errors in the three 
star-to-body angular measurements and in the range measurement (which includes 


I 


one of the star-to-body angular measurements and the Earth-Moon angular meas- 
urement) is' 

”1l/2 

rss = [(°Ax) 2 + (%^) 2 + ( a Az)j 

Actually, the quantities ox/or, dy/dr, and 'dz/'dr need not be calculated 
since the sum of the last terms in each of equations (C2) is approximately 
equal to the variance of the range error determined from equation (ll), where 
(for range to the Earth) 
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